Determination of the abundance of cosmic matter 
via the ceU count moments of the galaxy distribution 
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We demonstrate that accurate and precise cosmological information can be extracted from the 
cell count analysis of the 3D spatial clustering of galaxies once the second-order ratio between one- 
and two-point moments of the smoothed galaxy density distribution is analyzed. This probe does 
not require the calibration of any standard rod, the knowledge of the galaxy bias nor the modeling 
of galaxy redshift distortions. Using the spectroscopic Sloan Digital Sky Survey (SDSS) data release 
7 (DR7) galaxy sample, no cosmic microwave background (CMB) information, weak (flat) priors on 
the value of the curvature of the universe (fifc) and the constant value of the dark energy equation 
of state (w), we estimate the abundance of matter (Om) with a relative error of 11% (at 68% c.L). 
The method may be instrumental in searching for evidences of new physics beyond the standard 
model of cosmology and in planning future missions such as BigBOSS or EUCLID. 
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Determining the value of the constitutive parameters 
of the Friedman equations is a problem of great observa- 
tional and interpretative difficulty, but owing to its bear- 
ing on fundamental physics, is one of great theoretical 
significance. This challenge can be addressed through the 
development of performant cosmic probes, sophisticated 
error analysis schemes, and formidable observational pro- 
grams from both ground and space m. The Wilkinson 
microwave anisotropy probe (WMAP), for example, fix 
with impressive accuracy the parameters of a 'power-law 
ACDM' model [2|]. This is a cosmological model charac- 
terized by a flat geometry, a positive dark energy (DE) 
component fix with w — —1, and by primordial per- 
turbations that are scalar, Gaussian and adiabatic Q. 
Notwithstanding, other astrophysical probes are needed 
if we are to constrain deviations from this minimal model 
Different approaches, being sensitive to different sets 
of nuisance parameters and to different subsets of the full 
cosmological parameter set, provides consistency checks, 
lift parameter degeneracies, enable stronger constraints 
and, in the end, a safer theoretical interpretation. 

We characterize the inhomogeneous distribution of 
galaxies in terms of the local dimensionless density con- 
trast Sg{x) = pg{x)/{pg{x)) — 1, where Pg(x) is the co- 
moving density of galaxies in configuration space, and 
where (/9g(x)) denotes the selection function, the ex- 
pected mean number of galaxies at position x given 
the selection criteria of the survey. Since the galaxy 
distribution is a stochastic point process, a spheri- 
cal top-hat filter W of radius R is applied to gener- 
ate a continuous, coarse-grained observable (5g_ij;(x) — 
J (5(y)W^(|x - y|/i?)dy. The 1-point ^20,3,/?, = (5'^,i?(x))c 
and the 2-point Ku^gj^lr) = {Sg,fj{'x.)Sg,fj{'x. + r))c cumu- 
lant moments are the lowest-order non-zero connected 
moments of the probability distribution functional (PDF) 
of the field Sg^R{x) [l-Q. If the PDF is stationary and 



isotropic, the quantity Kii,g,i?(r)/K20,g,fl reduces to the 
ratio between the correlation function and the variance 
of the filtered field 



^g,R{r, p) 
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where we have explicitly emphasized the dependence of 
the observable on the set (p) of cosmological parameters. 
This comes trough because the LHS term of eq. ([T]) can 
be estimated from data only once a comoving distance- 
redshift conversion model is supplied. 

It is straightforward to predict the value of the second- 
order galaxy clustering ratio rjg j^. On sufficiently large 
cosmic scales R, matter fluctuations are small and 
described by the linear power spectrum A|^(A:,p) = 
47rylfc"=+'^r^(A;, p), where A is a normalization factor, Ug 
is the primordial spectral index, and is the transfer 
function [8] . Accordingly, the amplitudes of the second- 
order statistics for mass evolve as a function of redshift 
(z) and scale as ajf{z,p) = al{z — 0)D^{z)J^ii{p), and 
£_ji{r,z,p) = as{z — 0)^D^{z)Oji{r,p). Both these equa- 
tions are normalized at the scale rg = 8h~^Mpc, D{t) 
represents the linear growing mode Q, while the effects 
of flltcring are incorporated in the functions 



^k(p) 
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^^'^'^^ /+°°A|(fc,p)W^2(fc^g)rflnfc ^ ' 

where W is the Fourier transform of the window function. 

By analogy, we define the mass clustering ratio which 
is explicitly given by 

QRi-r, P) 
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As long as a linear matter power spectrum is assumed, rju 
is stable across cosmic time, i.e. it does not evolve as a 
function of redshift. Moreover, it is independent from the 
amplitude of the matter power spectrum normalization 
parameter A and it is also insensitive to the shape of the 
sample selection function, if the density gradient is nearly 
constant on scales ~ 2R. 

Interestingly, eq. ^ is also unaffected by redshift 
distortions, and, therefore, independent of their specific 
modeling, that is the ratio rif({r, p) between second-order 
statistics has identical amplitude in both configuration- 
and real-space. If the only net effect of peculiar velocity is 
the enhancement of the amplitudes of the density ripples 
in Fourier space, without any change in their phases or 
frequencies, then redshift distortions equally contributes 
to the numerator and denominator of eq. Q. For any 
given R and r, hence, the ratio r]fj{r, p) behaves as a 
cosmic 'standard of clustering'. 

The relationship between the mass and galaxy clus- 
tering ratios in configuration space follows immediately 
once we specify how well the overall matter distribution 
is traced by its luminous subcomponent on a given scale 
R. By adopting an arbitrary local non-linear biasing 

scheme, namely Sg,ji{x) = y]i3,o(^i/*')^fl(^)' ^^'^ 
limit (Tg j^ « 1, we obtain [Tf 

Vg,Rir,P) =r]R{r,p), (5) 

an identity that is independent of the specific value of the 
bias amplitudes bi. On scales R > 10/i~^Mpc and r — 
3R, the approximation ([5]) is good to better than 0.5%, 
and to better than 0.1% on scales larger than 15/i~^Mpc 
(for the same correlation length r = 3R). 

As for a classical Alcock and Paczynski test 0, the 
equivalence expressed by eq. ^ holds true if and only if 
the LHS and RHS are both estimated in the correct cos- 
mology. This is the principle that we exploit to constrain 
the value of key cosmological parameters. In this respect, 
we note that the RHS of eq. ([5]) is directly sensitive to the 
shape (ris) of the primordial power spectrum of density 
fluctuations, as well as to the parameters determining the 
size of the horizon scale at matter-radiation equality (the 
present-day extrapolated value of the matter (flmh'^) and 
baryon (fifj/i^) density parameters. The possible contri- 
bution of massive neutrinos is neglected in this analysis.) 
By contrast, the LHS term of eq. (0 probes the structure 
of the comoving distance-redshift relation, and, on top of 
fim, it is also sensitive to fix and w, i.e. to a wide range 
of DE models. Note that if distances are expressed in 
units of /i = i/o/lOO, the LHS is effectively independent 
from the value of the Hubble constant Hg. 

An additional merit of eq. ([5]) is the neat distinction in 
the physics that is brought about by its two terms. The 
LHS is a geometric probe since it is fully specified once 
the homogeneous expansion rate history of the universe is 
known. The nature of the RHS, instead, is fundamentally 
dynamical and it can be computed only when a specific 
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FIG. 1: Left panel: two-dimensional confidence limits on 
Qm and Q,a from a 'blind' analysis of the rCDM HVP sim- 
ulation This synthetic catalog of clusters of galaxies 
covers one octant of sky, extends over the redshift interval 
0.1 < z < 0.43, is comprised of ~ 10^ massive halos with an 
average space density of 6 x 10~*/i'^Mpc~^ and it was simu- 
lated using p' = (1, 0, —1, 21, 0, 1). Contours are plotted 
for C/£min < 2.3,6.17 corresponding to 68 and 95 per cent 
c.l. for a multivariate Gaussian distribution with 2 degrees 
of freedom. The relevant one-dimensional marginalized con- 
straint is shown in the inset. The clustering ratio %,fl(r, p) 
is estimated using R — 14/i~^Mpc and r — 3R. For the best 
fitting cosmology we measure r^g^n ir,pf) = 0.0294 ± 0.0017 
where the error is evaluated via 30 jackknife resampling of 
the data, excluding each time a sky area of 36 x 19.5 deg^. 
Dirac delta priors are taken of 0,bh^, Ho and n^, that are 
centered on the simulated values. Right panel: same as be- 
fore, but now contours represent the joint likelihood analysis 
of 40 independent ACDM mock catalogs simulating the SDSS 
LRG sample. Each of these simulations covers 120 x 45 deg^, 
extends over the redshift interval 0.16 < z < 0.43, is com- 
prised of ~ 60000 galaxies with an average space density of 
^ 8.9 X lO^'^^'^Mpc"'' and the input cosmological parameters 
are p' = (0.25, 0.75, -1, 70, 0.0196, 1). For the best fitting 
cosmology we find f?9,fl('', p^) ~ 0.0573 ± 0.0010, where the 
error, computed as the s. d. of the mean of 40 measurements, 
includes the contribution from cosmic variance. 

gravitational theory is provided. If the expansion rate of 
the universe is independently known, one can therefore 
test whether general relativity correctly describes the for- 
mation of correlated structures in the universe. Such a 
strategy complements other approaches exploiting red- 
shift distortion effects {e.g. [lO|). 

We estimate cr^ ^ and ^g.R and correct measurements 
for shot noise effects according to the method detailed 
in [3|- Here, we describe how we evaluate P(p|77(,^ij), 
the likelihood of the unknown set of parameters p = 
{Q.m,^x,w, Hi),Q,iih? ,ns) given the actual value of the 
observable ?7g,fl. Note that the analysis does not require 
the specification of further model parameters other than 
those quoted above. Since rjg ji is approximately Gaus- 
sian, the most likely set (p-^) is the one that minimizes 
the logarithmic posterior £ = — log P 

C = ^\oga^,i + — + (6) 
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FIG. 2: Upper panels: The galaxy clustering ratio 
rjg.ii{r,Qrn,^x ,w) (dotted line) reconstructed from SDSS 
DR7 data assuming R — 14/i~^Mpc and r = 37? is plotted 
as a function of Qa (left), f2m (center) and w (right), after 
fixing the remaining cosmological parameters, shown in the 
inset, to the values that minimize the logarithmic posterior 
((B)). Grey shading is the 1-sigma deviation around the mean 
and it is computed from 30 jackknife resampling of the data, 
excluding, each time, a sky area of 12 X 14 deg^ This un- 
certainty is in good agreement (< 10% difference) with the 
standard deviation displayed by the 40 SDSS-like simulations 
LasDamas. We also show the predicted scaling of the mass 
clustering ratio (solid line) in the corresponding cosmology. 
The correct set of cosmological parameters is the very one 
that makes the two curves converge to the same numerical 
value. For these data we find rig,R{r,pf) = 0.0502 ± 0.0036. 
Lower panels: two-dimensional marginalized constraints on a 
curved XCDM model in which both fix and w are allowed to 
vary. Gaussian priors are taken of Q^h^ — 0.0213 ±0.0010, of 
Ho = 73.8 ±2.4 and of = 0.96 ±0.014 from BBN [H, HST 
and WMAP7 2] determinations respectively. Contours 
are plotted for Lj Lmin < 2.3, 6.17. 



where — '^^^iVg-R ^ 'Hr)'^ ^ where N is the num- 
ber of estimates of cr^.i is the error on the ob- 
servable, TT describes any a-priori information about the 
PDF of p, and where B is a normahzation constant 
that can be fixed by requiring J Pdp — 1. Since eq. 
([5]) is free from look-back time effects, the analysis does 
not require slicing the sample in arbitrary redshift bins; 
only one estimation {N = 1) of eq. ^ is needed 
across the whole sample volume. We have recalculated 
the observable 77gjf(r, p) for each comoving distance- 
redshift model, i.e. on a grid {Qm,^x,w) of spacing 
[10~^, 5-10~^, 5-10~^]. As a consequence, the posterior C 
does not vary smoothly between different models because 
the number of galaxies counted in any given cell varies 
from model to model. However, since the computation 
of the observable rig^]i{3R,p) takes a limited amount of 
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FIG. 3: Two-dimensional marginalized constraints on Qrn and 
Ho obtained by fitting data with a curved XCDM model. 
Contours are plotted for £,/£rn,in < 2.3,6.17 Gaussian priors 
are taken of fi^/i^ = 0.0213 ± 0.0010 and Us = 0.96 ± 0.014 
from BBN ^ and WMAP7 [2] determinations respectively. 



time, shot-noise is the price we have decided to pay in 
order to obtain an unbiased likelihood hyper-surface. 

We test the performances of our method under realistic 
operating conditions via a 'blind' analysis of mock cata- 
logs that are characterized by widely different sets of ex- 
pansion rate and power spectrum parameters, mass trac- 
ers and radial selection functions. These are the rCDM 
Hubble Volume Project (HVP) simulation [llj and the 
LasDamas simulation [12]. An end-to-end analysis of 
clusters mock data (HVP) allows us to check for the pres- 
ence of insidious algorithmic biases that could arise when 
training a method to 'recognize' only a fiducial ACDM 
cosmology via a single class of biased tracer of the mat- 
ter clustering pattern, namely galaxies. The left panel of 
FIG. [T] shows that the input value Hm = 1 is statistically 
retrieved. By a joint analysis of 40 ACDM mock catalogs 
that incorporate all the observing selections of the SDSS 
luminous red galaxies (LRG) survey, i.e. the Oriana out- 
put of the LasDamas simulations, we have verified that 
the specific SDSS observing biases do not spoil our cos- 
mological inferences. The encouraging outcome of this 
analysis is presented in FIG. [U and shows that redshift 
data in a volume 40 times larger than that probed by the 
SDSS DR7 LRG sample (still - 10% of those that will be 
measured by surveys such as EUCLID J^] or BigBOSS 
[l6| ) can constrain fl^ precisely (~ 3.7%) and accurately 
(the input value is recovered). The strong sensitivity to 
the abundance of matter essentially arises because the 
zero order spherical Bessel function in eq. (3) filters in 
different portions of when is changed. If we bias 
high flm, the suppression of power in £,R{r, p) on a scale 
r is larger than that observed when estimating ^g,fl(r, p) 
using the corresponding wrong distance-redshift relation 
(see FIG. 2). Concerning DE parameters, and as for 

constraints are 



other clustering probes (e.g. 
not yet tight enough to yield a sensible measurement. 
Turning to real data, the geometry of the SDSS DR7 
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subsample 17[ is dictated by the need of sampling the 
galaxy distribution with cells of radius R as well as of 
correlating cell counts on scales r. The inferior limit 
on R is set by a theoretical constraint. Although not 
mandatory, the analysis is performed in the linear do- 
main, so as to avoid any phenomenological description 
of the matter power spectrum. The choice of R, is ad- 
ditionally conditioned by the practical requirement of 
minimizing the shot noise contribution in each cell, i.e. 
R > (47rp/3)~^^'^. In this work we adopt the scale 
R = 14/i~^Mpc. The inferior limit r = 3i? in the correla- 
tion scale, instead, is set to guarantee optimal accuracy 
in the approximation ([5]). At the opposite end, the larger 
the values of R and r, and the less informative are the 
inferred cosmological predictions. The sample that com- 
plies with these constraints covers the redshift interval 
0.15 < z < 0.43, a contiguous sky area of 120 x 45 deg^, 
is comprised of 62,652 LRG, and has a mean density of 
9.2 X 10-^h^Mpc~^. 

Without fixing neither the curvature of the universe 
(flat prior < \flk\ < 1) nor the quality of the DE 
component (flat prior —1/2 < w < —3/2), but tak- 
ing (strong) Gaussian priors of 0.bh? , Hq and Ug from 
big bang nucleosynthesis (BBN)ll8|, Hubble Space Tele- 
scope (HST)[i3], and WMAP [^respectively, data con- 
strain the overall matter abundance with a precision 
of 8% (see FIG. [H for the sake of completeness, by 
choosing a smoothing scale of 18/22/i~^Mpc we obtain 
0-270to;o27/0.255±0.038.) This figure improves by a fac- 
tor of two the estimate (fim — 0.259 ±0.039) obtained by 
analyzing, with the same priors, the large scale clustering 
(LSC) of the SDSS DR5 sample (which contains 25% less 
galaxies than DRT)^. It also improves by ^ 20% the 
constraint {flm = 0.24to;o24) of ^ obtained by combin- 
ing LSC results from the DR7 sample with the full like- 
lihood of the WMAP5 data [2l| and a strong HST prior 
[221] . The relative precision achieved with our method 
compares with that (^m = 0.294 ± 0.023) obtained by 
[i^ l by combining the LSC likelihood from the DR9 sam- 
ple (containing 5 times more LRGs than the DR7 cata- 
log) with the full likelihood of CMB data (WMAP7 0). 

Interestingly, the best fitting value of ftm still remains 
within the quoted 68% c.l. even when we relax some of 
the strong priors. With a (weak) flat prior on Us (in the 
interval [0.9,1.1]) we obtain Urn = 0.27lln;"3?, while if 
we additionally weaken the prior on flhh^ (flat in the in- 
terval [0,0.03]) we obtain = 0.255 ±0.040. Although 
the best fitting value of is weakly sensitive to changes 
in ilfjh^ and Ug, FIG. |3] shows that ftm degenerates with 
Hq when the HST prior is removed. If we relax also the 
strong prior on Hq (by assuming a flat prior in the in- 
terval [40, 100]km/s/Mpc), we find Qm = 0.22l°;Jt The 
stability of the best fitting central value fim is of even 
more interest if contrasted to CMB results showing that 
it is the combination flmh^ which is insensitive to the 
prior on the curvature of the universe. 



The relevant cosmological information contained in the 
SDSS DR7 data can be effectively retrieved by approxi- 
mating the galaxy clustering ratio rig^n in the likelihood 
expression Q via the fitting formula 

?7g,fl(3i?,p) = 0.00712-0.07882:p+0.0981x^-0.0538a;^ 

where R — 14/i^^Mpc, Xp = 1 + fp/8.48, and where 



fp — 'iwflx 



- Qk + 0.042 |wSlx [1-1(1 + w)w] + I 



-0.29 |u;r2x(14 ± 3w) - SwUxn^ 



in 

2 



} 



and by further assuming cr^ = 0.004, independently of 
the cosmological model. 

The advantage of probing cosmology via the ratio of 
second-order, one- and two-point moments of the dn dis- 
tribution is that only a minimum amount of cosmological 
hypotheses, and no astrophysical nuisance parameters at 
all, condition the results. Of great interest is the flexi- 
bility of the method, which can still be improved along 
several directions. Owing to its scale-free nature, the 
precision of the technique could be further improved by 
adopting a non-linear power spectrum in cqs. ^ and 
([3]), and by smoothing the over-density field on an even 
smaller scales R than that adopted throughout this anal- 
ysis. The only caveat being that mass fluctuations Sr 
in configuration and real spaces be simply proportional. 
Alternatively, the method could be applied also to large 
photometric redshift surveys. In this case, predictions 
of eq. dH) must be statistically corrected to account for 
the line-of-sight distortions introduced when estimating 
low-resolution distances from photometry. Finally, this 
probe enriches the arsenal of methods with which the 
next generation of redshift surveys such as EUCLID and 
BigBOSS will hunt for new physics by challenging all sec- 
tors of the cosmological model. 
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